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Abstract 

Two-particle angular correlations between trigger particles in the forward pseudorapidity 
range (2.5 < |rj| < 4.0) and associated particles in the central range (|rj| < 1.0) are mea¬ 
sured with the ALICE detector in p-Pb collisions at a nucleon-nucleon centre-of-mass en¬ 
ergy of 5.02 TeV. The trigger particles are reconstructed using the muon spectrometer, and 
the associated particles by the central barrel tracking detectors. In high-multiplicity events, 
the double-ridge structure, previously discovered in two-particle angular correlations at 
midrapidity, is found to persist to the pseudorapidity ranges studied in this Letter. The 
second-order Fourier coefficients for muons in high-multiplicity events are extracted after 
jet-like correlations from low-multiplicity events have been subtracted. The coefficients are 
found to have a similar transverse momentum (pj) dependence in p-going (p-Pb) and Pb- 
going (Pb-p) configurations, with the Pb-going coefficients larger by about 16 ± 6%, rather 
independent of p-\ within the uncertainties of the measurement. The data are compared with 
calculations using the AMPT model, which predicts a different pj and p dependence than 
observed in the data. The results are sensitive to the parent particle V 2 and composition of 
reconstructed muon tracks, where the contribution from heavy flavour decays is expected 
to dominate at pj > 2 GeV/c. 


© 2015 CERN for the benefit of the ALICE Collaboration. 

Reproduction of this article or parts of it is allowed as specified in the CC-BY-4.0 license. 


See Appendix |aJ for the list of collaboration members 






Forward-central two-particle correlations in p-Pb collisions 


ALICE Collaboration 


1 Introduction 

Measurements of correlations in A (p and At], where A (p and At] are the differences in azimuthal 
angle (c p) and pseudorapidity ( 77 ) between two particles, respectively, provide insight on the 
underlying mechanism of particle production in collisions of hadrons and nuclei at high energy. 

For such measurements in proton-proton (pp) collisions, jet production leads to a characteristic 
peak-like structure on the “near side” (at A (p « 0, A?] ca 0) and an elongated structure in A 77 
on the “away side” (at A (p ~ tt) OJ. In nucleus-nucleus (A-A) collisions, ridge-like structures 
extending over a long range along the A 77 axis emerge on the near and away sides, in addition to 
the jet-related correlations EHl4ll. The Fourier decomposition of the correlation in A (p at large 
Ar] is dominated by the second- and third-order harmonic coefficients V2 and V3, but significant 
harmonics have been measured up to vy, If 6 l 17.19l4l6ll. In A-A collisions, the v n coefficients are 
interpreted as the collective response of the created matter to the collision geometry and fluctua¬ 
tions in the initial state [IT71I181 . and are used to extract its transport properties in hydrodynamic 
models [H9l42T| . 

Long-range ridge structures on the near side (A (p « 0) were also observed in high-multiplicity pp 
collisions at a centre-of-mass energy y/s — 7 TeV [|22l and in proton-lead (p-Pb) collisions at a 
nucleon-nucleon centre-of-mass energy = 5.02 TeV Ii23l . Shortly after, measurements in 
which the contributions from jet fragmentation were suppressed by subtracting the correlations 
extracted from low-multiplicity events, revealed the presence of essentially the same long-range 
structures on the away side as on the near side in high-multiplicity events [24, i25j. Evidence 
of long-range double-ridge structures in high-multiplicity deuteron-gold (d-Au) collisions at 
y/^NN = 0.2TeV was also reported lf26Tl . By now, the existence of long-range correlations 
in p-Pb collisions is firmly established by measurements ll27j - [3Tll involving four, six or more 
particle correlations, with the lower-order correlations removed |[32l . demonstrating that the 
long-range ridges originate from genuine multi-particle correlations. Intriguingly, the transverse 
momentum dependence of the extracted v n [[271 (28, [30], and the particle-mass dependence of 
v n Q3I4351 are found to be qualitatively similar to those measured in A-A collisions. 

The similarity of the ridges in the pp, p-Pb, d-Au and A-A systems suggests the possibility 
of a common hydrodynamical origin 0614431 . However, whether hydrodynamical models can 
indeed be reliably applied to such small systems is under intense debate [|44l . Other proposed 
mechanisms involve initial-state effects, such as gluon saturation and extended color connec¬ 
tions forming along the longitudinal direction [14514491 or final-state parton-parton induced in¬ 
teractions ll50l - [54l . 

Further insight into the production mechanism of these long-range correlation structures may 
be gained by studying their 77 -dependence. A preliminary result Ii55l indicates a mild 77 depen¬ 
dence, but the measurement is limited to 1 77 1 <2. A similar magnitude of the two-particle cor¬ 
relation amplitudes in the Au-going and d-going directions at 2.8 < 77 1 <3.8 has also been re¬ 
ported in d-Au collisions at ^/snn = 0.2 TeV 11561 . Calculations for V 2 at large 77 (2.5 < 1 77 1 < 4) 
in p-Pb collisions at y/^NN = 5.02 TeV from a 3+1 dimensional, viscous hydrodynamical model 
and a multi-phase transport model (AMPT) predict a stronger 77 dependence, with about 50% 
and 30% larger +2 values on the lead nucleus side for the hydrodynamical and AMPT model, 
respectively Il57l . 

In this Letter, we report a measurement of angular correlations between trigger particles in the 
pseudorapidity range 2.5 < 1 77 1 < 4.0 and associated particles in the central range |rj| < 1.0 in 
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p-Pb collisions at y^VNN =5.02 TeV at the Large Hadron Collider (LHC). The trigger particles 
are inclusive muons, reconstructed using the ALICE muon spectrometer, and the associated par¬ 
ticles are charged particles, reconstructed by the ALICE central barrel tracking detectors. As in 
previous measurements [24, 33], the double ridge is extracted by subtracting the correlations ob¬ 
tained in low-multiplicity events from those in high-multiplicity events. Results for the second 
order Fourier coefficient for muons, {2PC, sub}, and the ratio of (2PC, sub} coefficient^] 
in the Pb-going (Pb-p) and p-going (p-Pb) directions are reported for high-multiplicity events, 
and compared to model predictions. The remainder of the Letter is structured as follows: We 
describe the experimental setup in Sec. [2], the event and track selection in Sec. [3} the analysis 
method in Sec. [4] and the evaluation of the systematic uncertainties in Sec. [5] Finally, in Sec. [ 6 ] 
we report the results, and compare them with model predictions. In Sec. [7] we conclude with a 
summary. 

2 Experimental setup 

In 2013, the LHC provided collisions between protons with a beam energy of 4 TeV and lead 
ions with a beam energy of 1.58 TeV per nucleon, resulting in a centre-of-mass energy of 
v /vnn = 5.02TeV. The beams were set up in two configurations: a period with the proton 
momentum in the direction of negative 77 in the ALICE coordinate system, denoted as p-Pb, 
followed by a period with reversed beams, denoted as Pb-p. Due to the asymmetric beam en¬ 
ergies, the nucleon-nucleon centre-of-mass reference system moves with a rapidity of 0.465 in 
the direction of the proton beam with respect to the ALICE laboratory system. Pseudorapidity, 
denoted by 77 , is given in the laboratory frame throughout this Letter. 

Details on ALICE and its subdetectors can be found in Refs. [|58ll59l . In the following, we give 
a brief summary of the components needed for the measurement reported in the Letter. 

Trigger tracks used in this analysis are detected in the muon spectrometer with an acceptance of 
—4.0 < 77 < —2.5. The muon spectrometer consists of a thick absorber of about ten interaction 
lengths (Ai), which filters muons in front of five tracking stations made of two planes of Cathode 
Pad Chambers each. The third station is placed inside a dipole magnet with a 3 Tm integrated 
field. The tracking apparatus is completed by a trigger system made of four layers of Resistive 
Plate Chambers placed behind a second absorber of 7.2 X\ thickness. This setup ensures that 
most of the hadrons in the acceptance are stopped in one of the absorber layers, providing a 
muon purity above 99% for the tracks used in this analysis. In p-Pb collisions, the trigger 
particle travels in the same direction as the p beam (p-going case), while in Pb-p collisions in 
the same direction as the Pb nucleus (Pb-going case). 

Associated particles in |rj| < 1.0 are reconstructed using the combined information from the 
Inner Tracking System (ITS) and the Time Projection Chamber (TPC), which are located inside 
the ALICE solenoid with a field of 0.5 T. The ITS consists of six layers of silicon detectors: 
two layers of Silicon Pixel Detector (SPD), surrounded by two layers of Silicon Drift Detector 
(SDD) and two layers of Silicon Strip Detector (SSD). SPD tracklets, short track segments 
reconstructed in the two SPD layers alone, are also used as associated particles. 

The V0 detector, consisting of two arrays with 32 scintillator tiles arranged in four rings each, 

1 Here, and in the following, “2PC” stands for “two-particle correlation” and “sub” for “subtraction”, and indi¬ 
cates the analysis technique with which the coefficients are measured. 
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Event (d/V ch /dT]) | ir) | <0 .5 
class pj > 0 GeY /c 
0=20% 35.8 ±0.8 

20-40% 23.2 ±0.5 

40-60% 15.8 ±0.4 

60-100% 6.8 ±0.2 

Table 1: VOS multiplicity classes as fractions of the analyzed event sample and the corresponding 
(d/V c h/dr]) ||fj|<o. 5 - The lAN^/Ar]) values are not corrected for trigger and vertex-reconstruction ineffi¬ 
ciencies, which are about 4% for non-single-diffractive events j6lil . mainly affecting the 80-100% lowest 
mulitiplicity events ll62l . Only systematic uncertainties are listed, since the statistical uncertainties are 
negligible. 

is used to generate the minimum-bias trigger and offline for multiplicity selection lf60ll . The 
detector covers the full azimuth within 2.8 < r\ < 5.1 (V0-A) and —3.7 < r\ < —1.7 (V0- 
C). The timing information of the V0 is also used for offline rejection of interactions of the 
beam with residual gas. In addition, two neutron Zero Degree Calorimeters (ZDCs) located at 
±112.5 m (ZNA) and —112.5 m (ZNC) from the interaction point are used in the offline event 
selection and as an alternative approach to define event-multiplicity classes. 

3 Event and track selection 

The online event selection used in this analysis is based on a combination of minimum-bias 
(MB) and muon trigger inputs. The MB selection uses the coincidence between hits in the V0- 
A and V0-C detectors and covers 99.2% of the non-single-diffractive cross section as described 
in H6TTI . Only approximately 5% of the MB events contain one or more tracks reconstructed in 
the muon spectrometer. In order to increase the number of recorded events, the presence of at 
least one muon above a pj threshold was required in addition to the MB trigger condition. Two 
different thresholds were used: a low-px threshold corresponding to about 0.5 GeV/c (/l-low- 
Pt) and a higher pj threshold corresponding to about 4.2 GeV/c (/i-high -pj). These thresholds 
are not sharp and the reported values correspond to a 50% trigger probability for a muon candi¬ 
date. The integrated luminosity collected with /i-high -pj triggers is 5.0 nb 1 in the p-Pb and 
5.8 nb -1 in the Pb-p periods. The /i-low -pj trigger class was downscaled by a factor 10-35 
depending on the data taking conditions, resulting in an integrated luminosity of 0.28 nb -1 in 
the p-Pb and 0.26 nb 1 in the Pb-p periods. 

The TPC and SDD detectors have longer deadtime compared to the muon spectrometer, the SPD 
and the V0. Therefore, they were read out only in a fraction of /i-low -pj events (about 25% in 
p-Pb and below 10% in Pb-p collisions). Both muon-track and muon-tracklet correlation results 
were measured in the p-Pb configuration. For Pb-p collisions, only muon-tracklet correlations 
could be studied due to the significantly lower number of triggers with the TPC in the readout. 

The primary-vertex position is determined using reconstructed clusters in the SPD detector as 
described in Ref. 11591 . Only events with a reconstructed vertex coordinate along the beam 
direction (z vtx ) within 7 cm from the nominal interaction point are selected. The probability of 
multiple interactions in the same bunch crossing (pileup) was dependent on the beam conditions 
and always below 3%. Pileup events are removed by rejecting triggers with more than one 
reconstructed vertex. 
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Fig. 1 : Parent particle composition of reconstructed muon tracks (left panel) and reconstruction effi¬ 
ciency for muons from pion and kaon decays relative to that for heavy flavor (HF) decay muons (right 
panel) from a detector simulation of the ALICE muon spectrometer. 


All events were characterized by their event activity, and sorted into event classes. As in pre¬ 
vious studies [ 24, 33], the event characterization was based on the signal in the V0 detectors. 
However, unlike before, both beam orientations were investigated in this Letter. Therefore, 
the signals from only two out the four rings of V0-A and V0-C detectors were combined to 
guarantee a more symmetric acceptance. On the V0-A side, the two outermost rings with an ac¬ 
ceptance of 2.8 < r] < 3.9, while on the V0-C side the two innermost rings with an acceptance 
of —3.7 < 77 < —2.7 were used. This combination is called VOS in the following. The definition 
of the event classes as fractions of the analyzed event sample and their corresponding average 
number of particles at midrapidity ((dAch/drj) r? |< 0 . 5 ), measured using tracklets as explained 
below, is given in Tab. |Tj 

Muon tracks are reconstructed in the geometrical acceptance of the muon spectrometer (—4 < 
T] < —2.5). The tracks are required to exit the front absorber at a radial distance from the 
beam axis, R a bs> in the range 17.6 < R a b s < 89.5 cm in order to avoid regions with large mate¬ 
rial density. The muon identification is performed by matching the tracks reconstructed in the 
tracking chambers with the corresponding track segments in the trigger chambers. Beam-gas 
tracks, which do not point to the interaction vertex, are removed by a selection on the product of 
the total momentum of a given track and its distance to the interaction vertex in the transverse 
plane. In the analysis, muons in the transverse momentum range 0.5 < pj < 4 GeV/c were 
considered. 

Reconstructed muons mainly originate from weak decays of n, K and mesons from heavy 
flavor (HF) decays. Because of the different pj distribution of the various sources and the ab¬ 
sorber in front of the spectrometer, which suppresses by design weak decays from light hadrons, 
the parent particle composition for the reconstructed muon tracks changes as a function of pj. 
The composition shown as a function of the reconstructed pj in the left panel of Fig. [T] was 
evaluated using full detector simulations based on the DPMJET Monte Carlo (MC) event gen¬ 
erator If63ll . The detector response was simulated using GEANT3 for particle transport [f64| . The 
composition of parent particles in the simulation differs by less than 10 % for the two beam con- 

2 Here, and in the following, pions and kaons refer to the sum of both charge states. Neutral particles are also 
considered in the case of kaons. 
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figurations. The reconstructed muons are dominated by light-hadron decays below 1.5 GeV/c, 
and by heavy flavor decays above 2 GeV/c. No significant multiplicity dependence was found. 
Similar conclusions are obtained using simulations with the AMPT generator Il65l . 

Without strong model assumptions, one cannot deduce the composition of parent particles from 
the measured muon distribution, and correct the data for muon decay and absorber effects. For 
comparison of the v ’2 data with calculations, however, only relative contributions of the parent 
species matter. In order to ease future model calculations, the reconstruction efficiencies for 
muons from pion and kaon decays relative to those for muons from heavy flavor decays are 
provided in the right panel of Fig. [j] as a function of the generated decay muon pj in different 
pseudorapidity intervals. Contributions from muon decays of other particles are significantly 
smaller than those for pions and can be ignored. The systematic uncertainty on the relative 
efficiencies was estimated to be less than 5%. 

Tracks reconstructed in the ITS and the TPC are selected in the fiducial region 1 77 1 < 1 and 
0.5 < pj < 4 GeV/c. The track selection used in this Letter is the same as in Ref. Il24l . 

Tracklet candidates are formed using information on the position of the primary vertex and the 
two hits on the SPD layers Il 66 ll . located at a distance of 3.9 and 7.6 cm from the detector centre. 
The differences of the azimuthal (A<Ph, bending plane) and polar (AO^, non-bending direction) 
angles of the hits with respect to the primary vertex are used to select particles, typically with 
p t > 50 MeV/c. Particles below 50 MeV/c are mostly absorbed by material. Compared to 
previous analyses lf 6 TI | 66 | a tighter cut in A(ph is applied (A<ph < 5 mrad) to select particles 
with larger pj and to minimize contributions of fake and secondary tracks to below 2.5%. 
The corresponding mean pj of selected particles, estimated from the DPMJET MC, is about 
0.75 GeV/c. 

4 Analysis 

The associated yield of tracks or tracklets per trigger particle in the muon spectrometer is mea¬ 
sured as a function of the difference in azimuthal angle (A (p) and pseudorapidity (At]). As in 
previous analyses 1124,1331, it is defined as 

Y _ 1 d^Aassoc _ S{Ar\,A(p) 

AtrigdArjdA^ B(Arj,A(p)’ 

in intervals of event multiplicity and trigger particle transverse momentum, p l T . The variable 
A tr ig denotes the total number of trigger particles in the event class and p\ interval, not corrected 
for single-muon efficiency. The signal distribution S(Ar] 1 A(p) — I/A lr igd 2 A sam e/dAr]dA(p is the 
associated yield per trigger particle for particle pairs from the same event, obtained in 1 cm-wide 
intervals of z vtx . A correction for pair acceptance and pair efficiency is obtained by dividing by 
the background distribution B(Ar],A(p) = a d 2 A m j xec j/dAr]dA(p. The background distribution 
is constructed by correlating trigger particles from one event with the associated particles from 
other events within the same event multiplicity class and 1 cm-wide z vtx intervals. The factor 
a is used to normalize the background distribution to unity in the A 77 region of maximal pair 
acceptance. The final per-trigger yield is obtained by calculating the average over the z vtx 
intervals weighted by A tr i g . 

In Fig. [2] the associated yield per trigger particle as a function of A (p and At] for muon-track 
correlations in p-Pb (left) and muon-tracklet correlations in p-Pb (middle) and Pb-p (right 
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Fig. 2: Associated yield per trigger particle as a function of Arj and A ip for muon-track correlations in 
p-Pb (left) and muon-tracklet correlations in p-Pb (middle) and Pb-p (right panels), measured in 60- 
100% (top row) and 0-20% (bottom row) event classes. The trigger particle (muon) range is 0.5 < p l T < 
1 GeV/c, the associated particle intervals are 0.5 < p\ <4.0 GeV/c for tracks and 0 < A<ph < 5 mrad 
for tracklets. Statistical uncertainties are not shown. 


panels), measured in 60-100% (top row) and 0-20% (bottom row) event classes is shown. 
In the low-multiplicity class (60-100%), the dominant feature is the recoil jet on the away 
side {%/2 < Acp < 3 k/2). While in previous two-particle correlation studies at midrapid¬ 
ity [24,'331 the away-side jet structure was mostly flat in At], from Arj = —1.5 to Arj — —5.0 
it decreases, as expected considering the ki nematics of dijets at large Arj. The near side 
(|A<p| < k/ 2) shows almost no structure in Acp and At], since it is sufficiently separated from the 
near-side jet peak at (A(p, Arj) — (0,0), so that no contribution from jets is expected. In the high- 
multiplicity (0-20%) class, the away-side jet structure is also visible, and the associated yields 
are considerably higher than for the low-multiplicity (60-100%) class. Moreover, in contrast to 
the low-multiplicity class, a near-side structure emerges, similar to that previously observed at 
lower pseudorapidities, revealing that the near-side ridge extends up to pseudorapidity ranges 
of 2.5 < \ri \ < 4. 

In order to isolate long-range correlations, we apply the same subtraction method as in previous 
measurements [24. 331 - Jet-associated yields have only a weak multiplicity dependence ll67l . 
thus the subtraction of the low-multiplicity event class removes most of the jet-like correlations. 
The per-trigger yield of the 60-100% event class is subtracted from that in the 0-20% event 
class, and the result is presented (labelled as y su b) in the top panels of Fig. [3j After subtraction, 
two similar ridges on the near and on the away side are clearly visible. 

The magnitude of the contributing long-range amplitudes is quantified by extracting the Fourier 
coefficients from the Acp projection of the per-trigger yield distribution, after the subtraction 
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Fig. 3: Top panels: Associated yield per trigger particle as a function of A<p and Arj for muon-track 
correlations in p-Pb (left) and muon-tracklet correlations in p-Pb (centre) and Pb-p (right) collisions 
for the 0-20% event class, where the corresponding correlation from the 60-100% event class has been 
subtracted. Statistical uncertainties are not shown. The trigger particle (muon) range is 0.5 < p\ < 
1 GeV/c, the associated particle intervals are 0.5 < p\ < 4.0 GeV/c for tracks and 0 < A<ph < 5 mrad 
for tracklets. Bottom panels: The same as above projected onto A(p. The lines indicate the fit to the data 
and the first harmonic contributions as explained in the text. 


of the low-multiplicity class, as shown in the lower panels of Fig. [3] In order to reduce the 
statistical fluctuations at the edges of the per-trigger yield distribution, the A (p projection is 
obtained from a first-order polynomial fit along Ar] for each A (p interval. In the p-Pb cases, the 
near- and away-side amplitudes are quite different, while in the Pb-p case the amplitudes on 
the near and away side are similar. The difference in the amplitudes of the near- and away-side 
ridge, which may be due to a residual jet contribution in the subtracted distribution, is taken into 
account in the systematic error evaluation, as explained in Sec. [5} 

The Fourier coefficients are then obtained by fitting T su b with 

ao + 2ai cos(A cp) + 2ci2Cos(2A(p) + 2«3COs(3A(p), (2) 

leading to ^ 2 /NDF values typically below 1.5. The relative modulation is given by V,,a{ 2PC, sub} 
pppyp, where b is the baseline of the low-multiplicity class (60-100%) estimated from the inte¬ 
gral of the per-trigger yield around the minimum. Assuming that the two-particle Fourier coef¬ 
ficient factorizes into a product of trigger and associate single-particle V 2 Il30l . the v„{2PC, sub} 
coefficients for particles reconstructed in the muon spectrometer are then obtained as 

v„{2PC, sub} = V,;a{2PC, sub}/y / Vj 7 c A {2PC, sub}, (3) 
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Assoc, tracks 

Assoc, tracklets 

Systematic effect 

p-Pb 

p-Pb 

Pb-p 

Ratio 

Acceptance (z vtx dependence) 

3-4% 

0-5% 

0-3% 

0-1% 

Remaining jet after subtraction 

4-10% 

5-14% 

1-2% 

3-15% 

Remaining ridge in low-multiplicity class 

1-4% 

1-6% 

0-2% 

2-8% 

Calculation of V 2 

0-1% 

0-1% 

1% 

0-2% 

Resolution correction 

1% 

0-1% 

0-1% 

0-2% 

Sum (added in quadrature) 

7-11% 

6-14% 

2-4% 

5-17% 


Table 2: Summary of main systematic uncertainties. The uncertainties usually depend on pj and vary 
within the given ranges. 

where V r ;? c A {2PC, sub} is measured by correlating only central barrel tracks (or tracklets) with 
each other (essentially repeating the analysis as in Ref. lf24lD . 

In this Letter, V2{2PC,sub} values for muons in the acceptance of the muon spectrometer are 
reported. Weak decays and scattering in the absorber of the muon spectrometer can cause 
the kinematics of reconstructed muons to deviate from those of their parent particles, and can 
influence the reconstructed V 2 , especially in case V 2 ,parent has a strong pj dependence. Since we 
cannot correct the measured V 2 for the species-dependent inefficiencies induced by the absorber, 
we denote the resulting coefficients by {2PC, sub} to indicate that the result holds for decay 
muons measured in the muon spectrometer. 

5 Systematic uncertainties 

The systematic uncertainty on {2PC,sub} was estimated by varying the analysis procedure 
as described in this section. The uncertainty on the ratio between the v^ (2PC, sub} in Pb-p and 
p-Pb collisions was obtained on the ratio itself, in order to properly treat the (anti-) correlated 
systematics between the p-Pb and Pb-p data samples. A summary is given in Tab.[2j 

The acceptance of the ALICE central barrel depends on the position of z vtx . To study its influ¬ 
ence on Vt (2PC, sub}, the analysis was repeated using only events with a reconstructed primary 
vertex within ±5 cm instead of ±7 cm from the nominal interaction point. The yield per trig¬ 
ger particle was not corrected for single track acceptance and efficiency of associated particles. 
Since V2(2PC,sub} is a relative quantity, it is not expected to depend on the normalization. 
This was verified in the case of the muon-track analysis, where good agreement was found be¬ 
tween the second-order Fourier coefficients obtained with and without single-track acceptance 
and efficiency corrections. Hence, no additional uncertainty was considered. 

As observed in previous analyses [24.132), the subtraction of the low-multiplicity class leads to 
a residual peak around (Ar /, A (p) « (0,0), possibly due to a bias of the event selection on the jet 
fragmentation in low-multiplicity events |[67l . The pseudorapidity gap [24}[25] used to calculate 
V„a was var ' c d from 1.2 to 1.0 and to 0.8 in order to estimate the contribution of the residual 
near-side short-range correlations. Due to the large gap in pseudorapidity between the ALICE 
central barrel and the muon spectrometer, this contribution does not affect the forward-central 
correlation. The effect of the bias introduced by the multiplicity selection was addressed on the 
away side by scaling the 60-100% multiplicity class. The scaling factor (/) is determined as 
the ratio between away-side yields in high- and low-multiplicity classes after the subtraction of 
the second-order Fourier component J67). This procedure was applied in the calculation of both 
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V nA and V£ A . The scaling factors were found to be larger in the case of p-Pb collisions (/ < 
1.40), compared to Pb-p (/ < 1.26), and tend to be lower for increasing pj. The difference 
with respect to the baseline results, for which no scaling (/ = 1) is applied, was taken as the 
systematic uncertainty. 

As previously reported ll67l . the contribution of the long-range correlations to the measured 
yields is not significant in low-multiplicity events. Still, their potential influence was addressed 
by changing the multiplicity range from 60-100% to 70-100% for the low-multiplicity class. 

To test the stability of the fit, the V 2 coefficient was calculated using a fit with only the first and 
the second Fourier components in Eq. |2j As another variation, the baseline b was calculated 
from a fit of the per-trigger yield in the low-multiplicity class using a Gaussian to model the 
shape of the aw ay-side ridge and a constant to estimate b. An equivalent approach, which 
makes use of the baseline of the high-multiplicity class B in V nA {2VC, sub} = a n /B , was also 
used, where B was estimated from the integral or from a parabolic fit of the correlation function 
around the minimum. Finally, the A (p projection was obtained from a weighted average instead 
of a first-order polynomial fit along At] for each A (p interval. 

The effect from the finite angular and momentum resolution of the muon spectrometer on 
i >2 (2PC, sub} was evaluated from a dedicated MC study with the measured V 2 as input distribu¬ 
tion, and resulted in a small correction of below 2%. The associated uncertainty was evaluated 
by varying the input V 2 by 50% at the lowest and highest measured points. 


6 Results 


The V2(2PC,sub} coefficients were measured for muon tracks in the p-going direction (p-Pb 
period) using both tracks and tracklets as associated central barrel particles, as described in 
Sec. [4 The (2PC, sub} coefficients obtained from the per-trigger yields of associated central 
barrel tracks agree well with those of associated tracklets, as shown in Fig. [4] as a function of 
muon pj. Since the two measurements probe different ranges in associated particle pj, the 
agreement is a consequence of trigger and associate V 2 factorization [|30l . In addition, good 
agreement was found between the (2PC, sub} obtained with different cuts on A(ph of associ¬ 
ated tracklets (inducing a change of average pj by about 20%). 


The p-going and Pb-going (2PC, sub} coefficients obtained using muon-tracklet correlations 
for the two different beam configurations (p-Pb and Pb-p) are reported in the left panel of 
Fig. 5 ] as a function of muon pj. The Pb-going V2(2PC,sub} (i.e. when the muon trigger 
particle travels in the same direction as the Pb nucleus) is observed to be larger than the p- 
going V 2 ( 2 PC,sub} over the measured pj range, but the two have a similar px-dcpcndcncc. 
To quantify the asymmetry, the Pb-going over p-going ratio for the (2PC, sub} coefficients 
is reported in the right panel of Fig. [5] as a function of muon pj. The ratio is found to be 
rather independent of pj given the statistical and systematic uncertainties of the measurement. 
A constant fit to the ratio adding statistical and systematic uncertainties in quadrature gives 
1.16 ±0.06 with a ^ 2 /NDF = 0.4. The analysis was also repeated using the energy deposited 
in the neutron ZDCs on the Pb-going side instead of the VOS amplitude for the event class 
definition. As discussed in detail in Il62ll . the correlation between forward energy measured in 
the ZDCs and particle density at central rapidities is weak in p-Pb collisions. Therefore, event 
classes defined as fixed fractions of the signal distribution in the ZDCs select different events, 
with different mean particle multiplicity at midrapidity, than the samples selected with the same 
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Fig. 4: Comparison of {2PC,sub} for —4 < I] < —2.5 extracted from nruon-track and nruon-tracklet 
correlations in p-Pb collisions at C^nn = 5.02 TeV. 



Fig. 5 : The v’2{2PC,sub} coefficients from muon-tracklet correlations in p-going and Pb-going direc¬ 
tions (left) and their ratio (right) for —4 < rj < —2.5 in p-Pb collisions at 0 vnn = 5.02 TeV. The data 
are compared to model calculations from AMPT. 


fractions in the V0 detector. Still, the {2PC, sub} values were measured to be similar, within 
25% of those extracted with VOS estimator. In addition, the asymmetry between Pb- and p- 
going v, {2PC. sub} was found to persist with similar shape and magnitude. The observed 
asymmetry may result from decorrelations of event planes at different rapidity Il68ll . 

The data in Fig. [5] can not be readily compared with existing predictions l(57Tl for a 3+1 di¬ 
mensional, viscous hydrodynamical model [[391 and the AMPT model with the string-melting 
mechanism enabled ll65ll . The model calculations were performed without taking into account 
the effect of the muon absorber, and represent the V 2 of primary particles, while as discussed 
in Sec. [3 the measured {2PC,sub} coefficients are reported for decay muons. Depending 
on particle composition and on the px-dependence of the parent particle V 2 distribution, the 
difference between primary particle V 2 and decay muon V 2 can be quite large. For example, at 
1 GeV/c, assuming the V 2 of the parent particles rises with pj like at mid-rapidity [f33ll . the mea¬ 
sured Vj (2PC, sub} for muons originating from decays of pions (kaons) would be ^20 (40)% 
larger than that of the parent pions (kaons). 
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Instead, in Fig. [5] we show a comparison of the data with AMPT model calculations performed 
with the same parameters as in If57ll . These calculations were performed at generator level, 
decaying primary particles into muons using the PYTHIA decayer ||69ll . The effects of the muon 
absorber were included by applying the pj and rj dependent relative efficiencies provided in the 
right panel of Fig. |T] Event characterization was done by mimicking the VOS criteria at particle 
level, i.e. by counting charged particles in 2.8 < q < 3.9 and —3.7 < q < —2.7. The V 2 values 
were obtained separately for muons decaying from pions, kaons and heavy-flavor hadrons, and 
otherwise performing the analysis in the same way as in data. We found the V 2 for HF muons 
to be consistent with zero within the generated statistics (5M events with a HF muon in the 
acceptance of the muon spectrometer for each period). Hence, for the inclusive V 2 , which is 
obtained by weighting the calculated V 2 with the relative yields in each decay channel, the V 2 for 
HF muons has been set to zero to reduce statistical fluctuations. In AMPT the factor / used to 
scale low-multiplicity class to eliminate the remaining jet contribution after subtraction, reaches 
values much larger than in the data, up to / = 2. Applying the scaling reduces the extracted V 2 
and consequently this choice constitutes the lower (upper) bound of the shaded area in Fig. [5] 
left (right), while the opposite bounds correspond to / = 1 (as used for the baseline result in the 
data). 

As shown in the left panel of Fig. [5J below pj < 1.5 GeV/c, where the inclusive muon yield 
is expected to be dominated by weak decays of pions and kaons, the calculation produces qual¬ 
itatively similar trends as observed in the data. However, quantitatively a different pj and i] 
dependence is found, visible in particular in the right panel of Fig. [5] At pj > 2 GeV/c, where 
the inclusive muon yield is dominated by heavy-flavor decays, the data may support a finite 
value for the V 2 of HF muons, or a drastically different composition of the parent distribution or 
their V 2 values in AMPT compared to data. Indeed, comparing predictions of AMPT to pion, 
kaon and D-meson yields measured at midrapidity 1:70. 71], muons from heavy-flavor decays 
would be underestimated by a factor 3-5 relative to pion and kaon decays assuming the same 
discrepancy between model and data at forward rapidity. A finite value for HF muon V 2 would 
be consistent with the emergence of radial flow in heavy-flavor meson spectra as predicted 
in [1721 . and has been recently measured in Pb-Pb collisions at ^nn — 2.76 TeV 11731 . 

7 Summary 

Two-particle angular correlations between trigger particles in the forward pseudorapidity range 
2.5 < 1 77 1 < 4.0 and associated particles in the central range |q | < 1.0 measured by ALICE are 
reported in p-Pb collisions at a nucleon-nucleon centre-of-mass energy of 5.02 TeV. The trigger 
particles are inclusive muons and the associated particles are charged particles, reconstructed 
by the muon spectrometer and central barrel tracking detectors, respectively. The composition 
of parent particles for the measured muons is expected to vary as a function of pj (Fig. |T|). A 
near-side ridge is observed in high-multiplicity events (Fig. [2]). After subtraction of jet-like cor¬ 
relations measured in low-multiplicity events, the double-ridge structure, previously discovered 
in two-particle angular correlations at midrapidity, is found to persist even in the pseudorapidity 
ranges studied here (Fig. [3]). The second-order Fourier coefficients for muon tracks are deter¬ 
mined assuming factorization of the Fourier coefficients at central and forward rapidity. The 
measurement in p-Pb collisions was performed in two different ways, using tracks or track- 
lets for particles at q < 1.0, yielding consistent results (Fig. [4]). The second-order Fourier 
coefficients for muons in high-multiplicity events were found to have a similar transverse mo- 
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mentum dependence in the p-going (p-Pb) and Pb-going (Pb-p) configurations, with the Pb- 
going coefficients larger by 16 ±6%, rather independent of pj within the uncertainties of the 
measurement (Fig. [5]). The results were compared with calculations using the AMPT model 
incorporating the effects of the muon absorber, showing a different pj and rj dependence than 
observed in the data. Above 2 GeV/c, the results are sensitive to the V 2 of heavy-flavor decay 
muons. Forthcoming model calculations should apply the relative efficiencies for muon decays 
from pion and kaons (provided in Fig. |T|) at generator level for detailed comparison with our 
data. Further measurements (e.g. of heavy-flavor muon yields or charged-particle V 2 at for¬ 
ward rapidity) will be needed to reduce the ambiguity between between muon parent particle 
composition and their V 2 - 
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